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ABSTRACT 

This is the second paper in a series studying the star formation rates, stellar com- 
ponents, metallicities, and star formation histories and evolution of a sample of blue 
compact galaxies. We analyzed spectral properties of 97 blue compact galaxies, ob- 
tained with the Beijing Astronomical Observatory (China) 2.16 m telescope, with spec- 
tral range 3580A - 7400A. We classify the spectra according to their emission lines: 
13 of the total 97 BCGs sample are non-emission line galaxies (non-ELGs); 10 have 
AG N- like emission (AGNs), and 74 of them are star- forming galaxies (SFGs). Emission 
line fluxes and equivalent widths, continuum fluxes, the 4000 A Balmer break index and 
equivalent widths of absorption lines are measured from the spectra. 

We investigate the emission line trends in the integrated spectra of the star-forming 
galaxies in our sample, and find that: 1) The equivalent widths of emission lines are 
correlated with the galaxy absolute blue magnitude Mb', lower luminosity systems tend 
to have larger equivalent widths. 2) The equivalent width ratio [Nii]6583/Hq is anti- 
correlated with equivalent width Ha; a relationship is given that can be used to remove 
the [Nil] contribution from blended Ha+[Nn]6548, 6583. 3) The [On], H/3, H7 and 
Ha fluxes are correlated; those can be used as star formation tracers in the blue. 4) 
The metallicity indices show trends with galaxy absolute magnitude and attenuation 
by dust; faint, low-mass BCGs have lower metallicity and color excess. 

Subject headings: galaxies: active - galaxies: compact - galaxies: general - galaxies: 
stellar content 



1. Introduction 

Blue compact galaxies (BCGs) are characterized by their very blue color, compact appearance, 
high gas content, strong nebular emission lines, and low chemical abundances (Kunth & Ostlin 
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2000; Ostlin et al. 2001). These properties are typical of unevolved systems, thus suggesting that 
BCGs should have suffered very few bursts of star formation during their lives and that some 
of them are probably experiencing their first burst. In a recent review, Kunth h, Ostlin (2000) 
argued that, despite a few remaining young galaxy candidates (like I Zw 18, SBS 0335-052), in 
most BCGs an old underlying stellar population does exist, revealing at least another burst of 
star formation (SF) prior to the present one (Papaderos etal. 1996; Kong & Cheng 2002). In 
addition, these properties make BCGs represent an extreme environment for star formation that 
differs from that in the Milky Way and in other quiescent nearby galaxies (Izotov & Thuan 1999; 
Izotov et al. 2001). Detailed studies of BCGs are important not only for understanding their 
intrinsic properties, but also for understanding of the chemical evolution of galaxies, for constraining 
models of stellar nucleosynthesis, for understanding star formation processes and galaxy evolution 
in different environments. 

To measure the current star formation rates, stellar components, metallicities, and star forma- 
tion histories and evolution of BCGs, we have prepared an atlas of optical spectra of the central 
regions of 97 blue compact galaxies in the first paper of this series (Kong & Cheng 2002; Paper I). 
Because we want to combine the optical spectra we obtained with Hi data to constrain simulta- 
neously the stellar and gas contents of BCGs, we selected most of our sample from Hi surveys by 
Gordon & Gottesman (1981). The spectra were obtained at the 2.16 m telescope at the XingLong 
Station of the Beijing Astronomical Observatory (BAO) in China. A 300 line mm -1 grating was 
used to achieve coverage in the wavelength region from 3580 to 7400 A with about 10A resolution. 

In the present paper, we provide measurements of emission line equivalent widths and fluxes, 
equivalent widths of absorption lines, 4000 A Balmer break index, as well as fluxes at several points 
of the continuum for our BCGs sample. We explore the trends in emission line fluxes and equivalent 
widths in the integrated spectra of SFGs in the sample. The absorption features and the continuum 
colors will be used to study the stellar population components and star formation history of blue 
compact galaxies. The emission line equivalent widths and fluxes will be used to determine the 
physical parameters of blue compact galaxies. 

The paper is organized as follows. In Section 2, we classify spectra according to their emission 
lines. Emission line equivalent widths and fluxes measurements are presented in Section 3. Section 
4 describes the continuum determination and the measurements of stellar absorption equivalent 
widths. In Section 5 we present an analysis of the emission line equivalent widths, line ratios and 
blue magnitudes of BCGs. Section 6 summarizes our conclusions. 

2. Spectral Classification 

Most of our BCGs sample was selected from the previous neutral hydrogen studies of blue 
compact galaxies by Gordon & Gottesman (1981). This study has focused on the Haro, Markarian, 
and Zwicky lists of galaxies and hence objects were selected on the basis of a blue color, UV- 
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excess or compactness, but not on the basis of emission line strength (Smoker et al. 2000). The 
optical spectral observations of these galaxies show a range in spectral properties; from galaxies 
with absorption line spectra to narrow emission line objects classified as SFGs. 

In order to study our sample galaxies in detail, we use emission lines to classify the sample 
spectra into three types: non-emission line galaxy (non-ELG), low-luminosity active galactic nuclei 
(AGN) and star- forming galaxy (SFG). Our classification scheme is outlined below. 

Because the Ha recombination line is easily detected in optical spectrum and only weakly 
affected by dust and underlying stellar absorption, we first separated the spectra into two broad 
categories, emission line galaxy and non-emission line galaxy, using the Ha recombination line. 
When Ha is detected in emission, we classified the galaxy as an emission line galaxy, otherwise 
as a non emission line galaxy. 13 of 97 BCGs spectra have no Ha recombination emission lines, 
and were classified as non-ELGs. Stellar Ha absorption lines are, in fact, prominent in all of these 
spectra. 

Next we classify the remaining 84 emission line spectra into active galactic nuclei and star- 
forming galaxies based on the next 3 steps: 

1. Because of their very large permitted line widths, Seyfert 1 galaxies are easily recognizable. 
We classify 5 galaxies (Mrk 335, Mrk 352, Mrk 6, I Zw 26, and Mrk 50) with broad emission 
lines in their spectra as Seyfert 1 galaxies. 

2. To discriminate the Seyfert 2 from star- forming galaxies, we consider the intensity ratio 
(corrected for reddening, see Sec. 3.4) R 23 = ([Oil] 3727+ [O m]4959, 5007) /H/3 (Dessauges 
- Zavadsky et al. 2000) to the remaining 78 galaxies except for VII Zw 631. Dessauges - 
Zavadsky et al. (2000) showed empirically that galaxies with log R23 ~ 1.1 are mostly Seyfert 
2. Based on this diagnostic, we classify 4 galaxies, Mrk 1, Mrk 622, Mrk 198, Mrk 270, as 
Seyfert 2. The classification of Seyfert 1 and Seyfert 2 agrees with the results of Veron-Cetty 
& Veron (2001). 

3. When neither [On] 3727 nor [O iii]4959,5007 are available, we use the flux ratio of Ha and 
[Nn]6583 to identify an AGN (Carter et al. 2001). If log([N n]6583/Ha) < -0.25, the galaxy 
is classified as a star-forming galaxy (SFG, Hn-like); otherwise it is classified as an AGN. 
This allows us to classify VII Zw 631, which has log([N n]6583/Ha)= -0.07, as an AGN. 

To summarize, our sample consists of 13 non-ELGs; 10 AGNs; and 74 SFGs. 
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Fig. 1. — Example of Gaussian fit for the blended lines Ha+[N n]6548, 6583. Three narrow Gaussian 
components are used. The crosses are observed spectral data. 
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3. 



Emission line measurements 



3.1. 



Measurements 



The prominent spectral features of SFGs and AGNs include some commonly found strong 
emission lines, such as [Oll]3727, H/34861, [O m]4959,5007, Ha6563, [Nn]6583, [Sll]6716, 6731, 
and some less common emission lines, such as H74340, [O m] 4363, Hei5876, and [Oi]6300. 

The rest-frame equivalent widths, EWs, and integrated fluxes, F, of the emission lines were 
measured by direct numerical integration, using the SPLOT program in IRAF. The continuum 
levels and integration limits for the lines were set interactively, with repeat measurements made in 
difficult case. For single emission lines such as H/34861, [Oni]5007, direct integral methods were 
used. This method allows the measurement of lines with asymmetric shapes (i.e. with deviations 
from Gaussian profiles). For blended lines such as Ha, [Nil] 6548,6583, and the [Sn]6716, 6731 
doublet, we used the Gaussian deblending program of SPLOT. In Figure 1, as an example, we 
show the three narrow Gaussian components to fit of Ha, [Nn]6548,6583 of III Zw 43. Note that in 
these blended cases, the lines are only partly blended. The interactive method allows us to control 
by eye the level of the continuum, taking into account defects that may be present around the line 
measured. It does not have the objectivity of automatic measurements, but it does allow us to 
obtain reliable, accurate measurements. 

The equivalent widths of various emission lines are listed in Table 1, for all SFGs and active 
galactic nuclei. The objects are ordered by increasing right ascension at the epoch 2000 (02000)- 
Column (1) lists the galaxy name (same as Table 1 of Paper I). Columns (2) - (9) list the equivalent 
widths of the commonly found emission lines. Columns (10) - (13) list the equivalent widths of 
less commonly found emission lines. The second line for each entry lists an estimate of the error 
(see in Sec 3.2). We use the convention that positive equivalent widths denote emission to conserve 
space and improve readability. A dash in the table indicates either that the corresponding segment 
of the spectrum is lacking or that the spectrum was too noisy in the region to give a reliable value 
of equivalent widths. We have chosen an equivalent width of 1.0 A as the lower limit for true 
detection. The observed emission line fluxes (the Galactic foreground reddening were corrected, 
see Kong & Cheng 2002) are listed in Table 2. 



For measurements of emission lines and absorption lines where the slope and curvature of 
the continuum are well defined, the main sources of random errors in the flux and equivalent 
width measurements are the uncertainty of the overall height of the continuum level, the individual 
intensity points within the interval of integration, the signal-to-noise ratio of the continuum, and the 
uncertainty in the choice of the best-fitting profile parameter. To estimate la standard deviations 
of emission lines, we followed the method outlined in Tresse et al. (1999), based on the formulae 



3.2. 



Standard deviations 
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Fig. 2. — Log of the detection level of equivalent widths, EW i.e. log (EW /a(EW)), versus 
log EW for 4 emission lines. The emission line name, the number (SFG + AGN, EW > 1.5A) 
of plotted data and the 2a, 3a detection levels (dashed lines) are indicated in each panel. The 
plotting symbols are coded according to spectral classification, the asterisks correspond to SFGs, 
the triangles to AGNs. 
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of propagation of errors and Poisson statistics. The derivation of error formulae can also be found 
in Longhetti et al. (1998). 

The error ap in the flux F of an emission line can be expressed as (Tresse et al. 1999): 

a F = a c D^2N~ + EW/D (1) 

The error gew with equivalent width EW of the line can be expressed as: 

EW I 

a EW = —a c D^fEW/D + 2N ptx + {EW/Df/N pix (2) 

where a c is the mean standard deviation per pixel of the continuum on each side of the line, D is 
the spectral dispersion in A per pixel (for our spectral sample, D=4.8) and N p i x is the number of 
pixels covered by the line. Because the signal-to-noise ratio for each pixel was estimated by scaling 
the continuum variance according to Poisson statistics, these are not exactly the formal statistical 
la errors of F or EW. Our approximation slightly overestimates the errors, but in our analysis, 
we are mainly interested in the consistency of the estimation of errors. 

In Figure 2, we plot a logarithm of the detection level of equivalent widths, log (EW/a(EW)), 
versus log EW for 4 commonly found emission lines, [On] 3727, H/3, [Om]5007, and Ha with 
equivalent widths above 1.5A. The result shows our equivalent widths limit is at a 3d confidence 
level for those commonly found emission lines, and the typical uncertainty in these measurements is 
less than 10%. For those less common lines, such as H74340, [O i]6300, the measurements typically 
have confidence levels > 2a, and the typical uncertainty in these measurements is about 20%. 



3.3. Stellar absorption correction 

Fluxes of emission lines will be used to determine the internal reddening of emission line regions, 
the star formation rate, and the element abundance of galaxies. It is known that the measurements 
are an underestimate of the real flux of the spectral lines, because of the underlying absorption 
component. To correct the underlying stellar absorption, some authors (such as Popescu & Hopp 
2000) adopt a constant equivalent width (1.5 — 2 A) for all the hydrogen absorption lines. Because 
the real value of the absorption equivalent width is uncertain and dependent on the age of star 
formation burst and star formation history (Izotov et al. 1994; Gonzalez Delgado et al. 1999), 
the other usual correction for the contamination by stellar absorption lines assumes absorption 
equivalent widths, and iterates until the color excesses derived from Ha/H/3, H/3/H7, and H/3/H5 
ratios converge to the same value (Izotov et al. 1994). 

To derive the absorption equivalent width for hydrogen lines, we have applied an empirical 
population synthesis method, which uses observed properties of star clusters as a base (Cid Fer- 
nandes et al. 2001), to our BCG spectra. This empirical population synthesis method can give the 
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synthetic stellar population spectrum, so we can measure these underlying stellar absorption fea- 
tures for hydrogen lines. A full description of this application and equivalent widths of underlying 
stellar absorption lines will be presented in a forthcoming paper. 

Emission line fluxes of Ha, H/3, and H7 are corrected for this underlying absorption effect as 
follows: 



HZ = C(l + EW^/EWS^), (3) 

where F{™ and F^ c s (see in Table 2, were corrected for Galactic extinction) are respectively 
the absorption corrected and the observed emission line fluxes, and EW{£* and EW^ (see in 
Table 1) are, respectively, the equivalent widths of the underlying stellar absorption line and of the 
observed emission line. 



3.4. Dust attenuation corrections 

The extinction of interstellar dust in SFGs modifies the spectra of these objects. It is necessary 
to correct all observed line fluxes for this internal reddening. The most widely method used to 
correct the emission line spectra for the presence of dust is based on the relative strengths of low 
order Balmer lines. In order to have an internally consistent sample, we applied this method to 
each of our objects, using only the ratio of the two strongest Balmer lines, Ha/H/3. 

We used the effective absorption curve t\ = Ty (A/5500A)- - 7 , which was introduced by Chariot 
& Fall (2000). The color excess arising from attenuation by dust in a galaxy, E l ^ v , can be written: 

Egt v = A v /R v = 1.086t v /R v , (4) 



= ln[F(H«)/F(H/3)] - ln[J(Ha)/J(H/?)] 
TV [(A Ha /5500)-°-7 - (A H/3 /5500)-°-7] 1 ' 

where I (Ra) / 1 (Hf3) is the intrinsic Balmer flux ratio, F (Ra) / F (R/3) is the observed Balmer flux 
ratio (were corrected for Galactic extinction and underlying stellar absorption), Ty is the effective 
F-band optical depth. X Ha = 6563A, A H/ 3 = 486lA, , and Ry = 3.1. We adopted I(Ra)/I(Rf3) 
= 2.87 for SFGs , and I(Ra)/I(Rf3) = 3.10 for the AGN-like objects (Dessauges - Zavadsky et 
al 2000). For 4 galaxies, I Zw 18, Haro 43, II Zw 70, and, I Zw 123, their observed flux ratio, 
F(Ha)/F(H/3), are less than the theoretical flux ratio 2.87, we set E l ^ v to be zero. The results of 
color excesses E % ^ v are listed in the last column of Table 3. 

The value of the color excess was then applied to the observed spectrum, and the final, intrinsic 
line fluxes relative to H/3 for each galaxy can be expressed as: 
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J ( A ) = F ( A ) p rv[(A/5500)- - 7 -(A H/3 /5500)- - 7 ] / fi x 

where F(X) and /(A) are the dust-obscured (observed) and intrinsic line fluxes, respectively. The 
attenuation-corrected line intensities relative to H/3 are given in Table 3 for each star forming galaxy 
and AGN. The objects are ordered by increasing right ascension at the epoch 2000 (02000)- Column 
(1) lists the galaxy name, columns (2) - (12) list the line intensities relative to H/3, column (13) 
list the intrinsic flux of H/3, which was corrected for both internal and Galactic extinction, and 
underlying stellar absorption. Column (14) lists the color excesses E l ^ v for each galaxy. 



3.5. Comparison with previous studies 

Seven galaxies in our BCG sample — III Zw 43 (0211+038), II Zw 40 (0553+033), Mrk 5 
(0635+756), I Zw 18 (0930+554), Haro 4 (1102+294), Haro 29 (1223+487), I Zw 123 (1535+554) 
- have been observed previously by Izotov, Thuan & Lipovetsky (1997, ITL97), Izotov & Thuan 
(1998, IT98), and Guseva, Izotov, & Thuan (2000, GIT00), with the Ritchey-Chretien spectrograph 
at the Kitt Peak National Observatory (KPNO) 4 m telescope, and with the GoldCam spectrograph 
at the 2.1 m KPNO telescope. These high signal-to-noise ratio spectrophotometric observations 
allow us to test the quality of our data. We perform a detailed comparison of these previous studies 
in this subsection. 

In III Zw 43, GIT00 did not detect [O ill] 4363 line, the [S 11] 6731 line intensity ratio is about 
10% higher, and the other emission line ratios are in good agreement with ours. For II Zw 40 in 
GIT00, Haro 29 in ITL97, our data are in fairly good agreement with these studies. In Mrk 5, 
our [Oi]6300 and [Sn]6717 line intensities are about 25% higher, and [On] is ~ 13% lower than 
that in IT98. In I Zw 18, our H7, [Oiii]4363, [Sn] line intensities are stronger, but Hei5876 is 
weaker. In Haro 4, some less strong lines are not good agreement with IT98. Finally for I Zw 123, 
the agreement is not good as the other galaxies, our I(A)//(H/3) data have large differences with 
ITL97, but the F(X)/F(R(3) are in good agreement with ITL97. 

We now display this comparison in a more visible form in Figure 3. The horizontal axis 
represent different spectral lines, the vertical axis shows the differences between our line intensities 
([/(A)//(H/3)]our ) and the values of GIT00, IT98 and ITL97 ([/(A)//(H/3)]i T ). We found, our 
line intensity ratios are in good agreement with these previous studies for most spectral lines of 
most galaxies, the difference between our sample and these studies is less than 10% for those strong 
emission lines, and less than 15% for those less strong lines, such as H74340JO iii]4363,[0 i]6300 of 
most galaxies. 

The observed fluxes of H/3 in our data are larger than those in previous studies, the explanation 
could be: 1) the data in these previous studies were not corrected for the Galactic extinction; 2) 
Our slit width is larger than that of previous studies; 3) the position angle of slit is different between 
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ours and those previous studies. We will discuss the slit effect and derive an aperture correction 
for each galaxy in a future paper. 

4. Continuum and absorption line measurements 

The equivalent widths of absorption lines and the continuum colors provide information about 
the stellar populations and physical parameters of galaxies. One of our goals is to study the star 
formation history and chemical evolution of BCGs. To this end, we also determined a pseudo- 
continuum at selected pivot points and measured the equivalent widths of 7 absorption features, 
integrating the flux within each window between the pseudocontinuum and the spectrum. The 
absorption features and continuum points are chosen based on the population synthesis method 
(Cid Fernandes et al. 2001) that will be used in a forthcoming paper. 

4.1. Continuum measurements 

In order to derive a continuum, we have measured the flux values at 9 pivot points, 3660, 4020, 
4510, 4630, 5313, 5500, 6080, 6630, 7043 A which were chosen to avoid regions of strong emission or 
absorption features (Bica 1988, Kong & Cheng 1999, Saraiva et al. 2001). The corresponding fluxes 
were determined as averages in 20 A bins centered in the listed wavelengths. The determination 
of the continuum has been checked interactively, taking into account the flux level, the noise and 
minor wavelength calibration uncertainties as well as anomalies due to the presence of emission 
lines. The excellent quality of the spectra allowed a precise determination of the continuum in the 
majority of cases. 

In addition, 3 point flux values (3784, 3814, 3918 A) were measured, which were necessary for 
the determination of the continuum in galaxies with strong contribution of late B to F stars which 
present several high-order Balmer absorption lines in this region (Bica et al. 1994). Due to the 
crowding of the absorption lines, it is difficult to make automatic measurements. We thus selected 
the highest value of A3784, A3814 and A3918 A fluxes to represent the continuum in these spectral 
regions. 

Figure 4 illustrates the application of the method to two of the sample spectra. The spectrum 
in Fig. 4a corresponds to the nucleus of the non-emission line galaxy II Zw 82. The spectrum is 
purely stellar, so it is straightforward to measure the pseudocontinuum. Fig. 4b shows the nuclear 
spectrum of Haro 1, a star-forming galaxy. The continuum is also well defined. Overall, we find 
that the method works well for most spectra, with the exception of the nuclear regions of AGNs 
(see end of this Section). 

The specific continuum wavelengths and corresponding fluxes for the non-ELGs and SFGs 
are shown in Table 4. The second line for each galaxy entry lists the errors with continuum 
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Fig. 3. — Difference between our line intensity ratios [J(A)/J(H/3)]oUR and these of previous studies 
[I (\) / 1 (H(3)]it . The plotting symbols represent different galaxies, the dash lines outline the 15 per 
cent error window. 
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Fig. 4. — Illustration of the continuum determination procedure for the nuclear spectra of a) the 
non-emission line galaxy II Zw 82, and b) the star-forming galaxy Haro 1. The continuum pivot 
points are marked by the filled triangles. Vertical lines at the bottom of each panels indicate the 
wavelength windows used to measure the equivalent widths of Ca K and H, H<5, CN, G band, Mg 
and Na lines (See in Table 5). 
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measurements, calculated as the rms deviation from the average continuum flux. 

4.2. Absorption line measurements 

The spectrum of a galaxy is produced by the sum of the spectral characteristics of its stellar 
content (Weiss et al. 1995). Observations of the integrated spectra of galaxies can be used to 
determine the distribution in age and metal abundance of the stellar population in these systems 
and hence to determine their epoch of formation and subsequent star formation history (Arimoto 
& Yoshii 1986). To this aim, some strong, easily identifiable absorption features in our observed 
spectral range are measured, which include some age and metallicity sensitive absorption features. 
The absorption line names and adopted spectral wavelength windows are shown in Table 5. 

The rest-frame equivalent widths of these absorption features were automatically computed 
by summing the observed fluxes below the continuum level, which itself is estimated by fitting a 
straight line to the fluxes in the above continuum regions. The equivalent widths of the absorption 
features for these non-ELGs and SFGs are presented in Table 6. The second line for each galaxy 
entry lists the errors on equivalent widths, which are computed from equation (2) in this paper. 
Figure 5 shows that logarithm of the detection level of EW versus log EW for 4 absorption lines, 
such as Ca K, Ca H, R5 and Mgl+MgH when its EW > 1.5A. Our absorption feature EW limit 
is at a 3a confidence level for the absorption lines of most galaxies. The mean uncertainty in these 
measurements is about 10% for those non-ELGs and about 15% for those SFGs. 

4.3. 4000 A Balmer break index measurements 

As well as the absorption features, we also measured the 4000 A Balmer break. It is the 
strongest discontinuity in the optical spectrum of a galaxy and arises because of the accumulation 

Table 5: Wavelength windows used to measure the equivalent widths of the absorption lines. 
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Fig. 5. — Log of the detection level of equivalent widths, i.e. log [EW/a(EW)], versus log EW for 4 
absorption lines. The absorption line name, the number (SFG + non-ELG, EW > 1.5 A) of plotted 
data and the 2a, 3a detection levels (dashed lines) are indicated in each panels. The plotting 
symbols are coded according to spectral classification, the asterisks correspond to star-forming 
galaxies (SFGs), the triangles to non-ELGs. 
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of a large number of spectral lines in a narrow wavelength region (Bruzual 1983). The main 
contribution to the opacity comes from ionized metals. In hot stars, the elements are multiple 
ionized and the opacity decreases, so the 4000 A break will be small for young stellar populations 
and large for old, metal-rich galaxies (Kauffmann et al. 2002). 

We use the definition using narrower continuum bands than Bruzual's, which was introduced 
by Balogh et al (1999). The principal advantage of the narrow definition is that the index is con- 
siderably less sensitive to reddening effects. It is defined as the ratio of the average fluxes (for 
frequency unit) measured in the spectral ranges 4000-4100A and 3850-3950A: D4000vn=F„ [4000- 
4100A]/TV[3850-3950A]. The D4000vn index is simply a flux ratio and, hence the error is deter- 
mined from standard propagation techniques. The D4000vn index value and its error for these 
non-ELGs and SFGs are presented in the last column of Table 6. 

While for non-ELGs and SFGs the placement of the continuum is straightforward, this is not 
the case in the nuclear regions of most AGNs, where the numerous broad lines and intense non-stellar 
continuum complicate the analysis. The continuum points are impossible to determine accurately, 
and the equivalent widths of the absorption lines cannot be measured accurately. Therefore, for 10 
AGNs, we only measured the integrated fluxes, F, and rest-frame equivalent widths, EWs, of the 
emission lines. We did not measure the continuum fluxes and equivalent widths of the absorption 
lines for those AGNs. 

5. Emission line analysis for SFGs 

The main goal of this series of papers is to measure the current star formation rates, stellar 
components, metallicities, and star formation histories and evolution of BCGs. For this purpose, 
we are mostly interested in the SFGs. Therefore, in present section, we do not consider the galaxies 
with Seyfert nuclei and the non-emission line galaxies in our BCG sample. In the following, the 
sample will refer to the 74 SFGs. We investigate the emission line trends in spectra of those SFGs 
in this section. 



5.1. Equivalent widths of emission lines 

5.1.1. EWs versus Mb 

It is interesting to explore how the equivalent widths of emission lines depend on the galaxy 
absolute blue magnitude Mb (Column 6 of Table 1 in Paper I). The equivalent widths of [O ll]3727, 
H/3, and Ha are well correlated with Mb, the other emission lines are also correlated with Mb, but 
the spread in equivalent widths at a given luminosity is large. Lower luminosity systems tend to 
have larger equivalent widths for most of emission lines, except for [Nll]6583. 
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In the top panel of Figure 6, we plot the Ha emission equivalent width as a function of Mg for 
the SFGs in our sample. A pronounced trend towards larger equivalent widths at lower luminosities 
can be found, galaxies with the strongest Ha lines are of low luminosity. EW(Ha) is the ratio of 
the flux originating from UV photoionization photons (< 912 A) to the flux from the old stellar 
population emitted in the rest-frame R passband, which forms the continuum at Ha. Thus, a 
large equivalent widths is due either to a large UV flux (or B absolute magnitude since they are 
correlated), or to a low continuum from old stars. In either case, this implies a blue continuum 
color. Hence, the observed trend of larger EW(Ha) for fainter galaxies implies that the faint SFG 
population is dominated by blue galaxies, while the bright SFG population is dominated by redder 
galaxies. 

In the bottom panel of Figure 6, we plot the [Nil] 6583 emission line equivalent width as a 
function of Mb for the SFGs in our sample. Its equivalent width behaves in the opposite way, lower 
luminosity systems tend to have smaller equivalent widths. Such a trend has also been found in 
other studies of nearby galaxies ( Jansen et al. 2000) . The global behavior of [N n] 6583 EW reflects 
intrinsic differences in the nitrogen abundance in BCGs, on average luminous BCGs are likely to 
be enhanced in nitrogen abundance. This suggests that, in faint, low-mass, BCGs, nitrogen is a 
primary element, whereas in brighter, more massive BCGs it comes from a secondary source. 

5.1.2. EW([Nll])/EW(Ha) versus EW (Ha + [Nil]) 

In deep large optical surveys, low-resolution spectroscopy or narrowband Ha imaging is often 
used. Ha and [Nn]6548, 6583 lines are often blended, so it is important to recover the flux solely 
in Ha to measure for instance the Ha luminosity function, hence to derive a star formation rate. 
Figure 7a shows that the [Nn]6583/Ha EW ratio decreases as a function of EW(Ra). All the 
spectra in this figure have £W([Nll]) and EW(Ra) > 10 A, which can be measured very accurately. 
The [Nn]6583/Ha equivalent widths ratio is strongly correlated with EW(Ha). A least-squares fit 
of this relation yields: log EW ([N n]6583) / EW (Ra) = (1.01 ± 0.14) - (0.85 ± 0.07) log EW(Ra). 

Since EW( [Nil] 6583) = 3£W( [Nil] 6548), we also plot the relation 1.33 £W([N ll]6583)/ 
EW(B.a) versus £W(Ha) + 1.33£W([N n]6583) in Figure 7b. The trend is similar to that in Fig. 7a, 
log 1.33£W([N n]6583)/£W(Ha) = (1.36±0.20)-(0.91±0.09) log (W(Ha)+1.33£7W([Nn]6583)). 
Thus we can predict which value is expected for the ratio when observing the blend Ha + [N ll]6548, 
6583. For instance, if this latter, EW(Ha) + 1.33EW([N n]6583), is ~ 100 A, the ratio 1.33 
£W([Nll]6583) / EW(Ra) should be ~ 0.35. The value of the ratio [Nil] 6548,6583 / Ha, as 
determined by Kennicutt (1992), is usually taken to be 0.5 to remove the contribution of [Nil] to 
(Ha+ [Nil]) blended emission. This is slightly larger than the typical value for our star-forming 
galaxy sample trend, presumably because Kennicutt's sample contains a large fraction of early-type 
galaxies, which have systematically higher ratios (Tresse et al. 1999). 
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Fig. 6. — The logarithm of integrated Ha and [N n]6583 emission line EWs plotted versus absolute 
B filter magnitude. A tendency towards larger Ha EWs and smaller [N n] EWs in lower luminosity 
systems is seen. 
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5.2. Fluxes of emission lines 

Emission line fluxes are primary traces of the star formation rate (SFR) in galaxies (Kennicutt 
1998). For star- forming galaxies, the Balmer emission line luminosities scale directly with the 
ionizing fluxes of the embedded young stars, and this makes it possible to use the Balmer lines to 
derive SFRs in galaxies. Ha is the best line for such applications, but beyond z ~ 0.2 — 0.3, this 
line is redshifted into the near infrared. To find other tracers of SFR we analyze the fluxes of other 
emission lines as a function of the intrinsic Ha flux. 

Since the various prominent emission lines correlate with each other, any of them is likely to 
be a first order ranking indicator of SFR of star- forming galaxies, but the strongest correlations are 
found between [Oll]3727, H74340, H/3 and Ha. [On]3727 is the most useful star formation tracer 
in the blue. In Figure 8a, we show the flux of [On] 3727 as a function of the intrinsic Ha flux. We 
found, as expected, these two lines have a strong correlation. 

From purely astrophysical considerations, the most reliable star formation tracers in the blue 
should be the higher order Balmer lines, since the fluxes of these lines scale directly with the massive 
star formation and are nearly independent of the temperature and ionization level of the emitting 
gas. Figure 8b shows the relation between the fluxes of Ha and H7. A strong, roughly linear 
correlation between Ha and H7 is apparent. This correlation confirms that the H7 line can serve 
as a reliable star formation tracer in strong emission line galaxies, such as the SFGs in our BCG 
sample. In addition, the correlation between the intrinsic fluxes of Ha and H/3 is stronger than 
that between the fluxes of Ha and H7, H/3 is another good star formation tracer for SFGs. 

Recently, Chariot Sz Longhetti (2001) quantified the uncertainties in these SFR estimators. 
They found SFR estimates based purely on one of emission line luminosity of galaxies could be in 
error by more than an order of magnitude. On the other hand, with the help of other emission lines, 
these errors can be substantially reduced. Based on our high quality spectrophotometric data, we 
can derive star formation rate for each galaxy, using these different star formation rate estimators. 

5.3. Line ratios and metallicity 

Numerous studies have claimed the existence of a metallicity-luminosity relation in a vari- 
ety of classes of galaxies: dynamically hot galaxies, i.e. ellipticals, bulges, and dwarf spheroidals, 
dwarf Hn galaxies, irregular galaxies and spirals (Lequeux et al. 1979; Skillman, Kennicutt, h 
Hodge 1989; Stasihska & Sodre 2001). In Figure 9a-9c, we show the various metallicity indices, 
i? 23 =([Oii]3727+[Oiii]4959, 5007) /H/3 (Pagel et al. 1979), [Nn]6583/Ha (van Zee et al. 1998) and 
[Nil] 6583/ [On] 3727 (Dopita et al. 2000) as a function of the total absolute blue magnitude M B . 
The line flux ratios were corrected for both internal (using the value of E^^y) and Galactic extinc- 
tion, and underlying stellar absorption. Since the reddening correction becomes more uncertain for 
galaxies with small EW(Ra) and EW(R(3), we only use those objects that have EW(R(3) > 5 A 
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and thus the most reliable reddening corrections. We do find a good correlation between Mb and 
the metallicity indices [Nn]6583/Ha and R 23 . The correlation of the [Nil] 6583/ [Oil] 3727 index 
with Mb is statistically significant, but at a rather low level. These metallicity indices show clear 
trends with galaxy absolute magnitude, confirming that indeed there is a relation in blue compact 
galaxies between the overall metallicity of the star forming regions and the galaxy luminosity. The 
higher the galaxy absolute magnitude, the higher the heavy element content. This relation suggests 
that the metallicity of faint, low mass BCGs is low. 

We now examine whether there is a relation between the color excesses due to internal extinc- 
tion, E^y, and the overall metallicity of the galaxies. So far, there have been contradictory claims 
in this respect. Zaritsky et al. (1994) found no evidence for a systematic dependence between red- 
dening and abundance in a sample of 39 disk galaxies. In other contexts, Stasihska & Sodre(2001) 
found that the nebular extinction as derived from the Balmer decrement strongly correlates with 
the effective metallicity of the emission line regions of spiral galaxies. 

Figure 9d shows E l ^ v as a function of the metallicity indicator [N n]6583/Ha, which is less 
affected by the reddening correction. We find there is a clear correlation, log(I[Nii]6583/-^Ha)cor = 
-(1.06 ± 0.07) + (0.92 ± 0.17)£j^ y . [Nn]6583/Ha tends to be larger for larger values of E%i v , 
when E^iy > 0.1. Internal extinction indeed correlates with the overall metallicity of BCGs, 
especially among the galaxies with large E^y. Since the metallicity indices correlate with Mb, 
the correction between /[NH]6583/^Ha and E l gt v also suggest the internal extinction of brighter, 
more massive BCGs is higher. 



6. Summary 

In this paper, we have measured the fluxes and equivalent widths of emission lines, the fluxes 
at several points of the continuum, the 4000 A Balmer break index, and the equivalent widths of 
absorption lines for our BCGs sample. Then we have analyzed the fluxes and equivalent widths of 
emission lines for the star-forming galaxy subsample. 

Our main results are: 

1. We have measured the fluxes and equivalent widths of emission lines for BCGs, and estimated 
their errors. The typical uncertainties of the measurements are less than 10% for the brightest 
Ha6563, H/34861, [Oll]3727, [O m]4959,5007, [Nn]6583, [Sll]6717,6731 lines and less than 
20% for the weaker H 7 4340, Hei5876, [Oi]6300 and [Om]4363 lines. Our line intensity 
ratios are in good agreement with those of previous studies. The color excesses due to internal 
extinction was calculated based on the flux ratio of Ha/H/3. 

2. The equivalent widths of absorption features, the continuum colors and the 4000 A Balmer 
break indices of BCGs were measured. The mean certainty in these measurements is better 
than 15%. 
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3. The equivalent width of Ha emission line is correlated with Mb] lower luminosity systems 
tend to have larger equivalent widths. This can be explained by assuming blue galaxies dom- 
inate the faint SFG population, while redder galaxies dominate the bright SFG population. 

4. On average, luminous BCGs are likely to be enhanced in nitrogen abundance. This suggests 
that in faint, low-mass, BCGs nitrogen is a primary element, whereas in brighter, more 
massive BCGs nitrogen comes from a secondary source. 

5. A correlation is found between most of emission lines, the strongest correlations are found 
between [On] 3727, H7, H/3 and Ha. Besides [On] 3727, H7 and H/3 lines can serve as a 
reliable star formation tracers for strong emission line galaxies. 

6. Metallicity indices show clear trends with galaxy absolute magnitude, confirming that there 
is a relation in galaxies between the overall metallicity of the star forming regions and the 
galaxy luminosity. Faint, low-mass BCGs have lower metallicity and internal extinction. 
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